Submitted to ApJ 2006 January 3, Accepted 2006 March 20 

Preprint typeset using 1^1^^ style cniulatcapj v. 6/22/04 



o 

o 

(N 
>^ 



(N 



(N 
> 

(N 

O 

o 
6 



% 



THE SURVEY FOR IONIZATION IN NEUTRAL GAS GALAXIES- II. THE STAR FORMATION RATE 

DENSITY OF THE LOCAL UNIVERSE 

D.J. Hanish\ G.R. Meurer\ H.C. Ferguson-, M.A. Zwaan^, T.M. Heckman\ L. Staveley- Smith*, J. 

Bland-Hawthorn"', V.A. Kilborn*'", B.S. KoribalskL*, M.E. Putman^, E.V. Ryan-Weber**, M.S. Oey", R.C. 

Kennicutt, Jr.", P.M. Knezek"', M.J. Meyer^ R.C. Smith", R.L. Webster^^ M.A. Dopita^-\ M.T. Doyle", M.J. 

Drinkwater", K.C. Freeman'^^, J.K. Werk'^ 

Submitted to ApJ 2006 January 3, Accepted 2006 March 20 

ABSTRACT 

We derive observed Ho and R band luminosity densities of an Hi-selected sample of nearby galaxies 
using the SINGG sample to be l'^^ = (9.4 ± 1.8) x lO^^ ^170 erg s'^ Mpc'^ for Ha and I'j^ = (4.4 ± 
0.7) X 10^^ /170 erg s~^ A^^ Mpc~'^ in the R band. This R band luminosity density is approximately 
70% of that found by the Sloan Digital Sky Survey. This leads to a local star formation rate density 
of log(/9sFR [-^Q yr~^ Mpc~^]) = —1.80 f^'l^{random) ±0.03{systematic) + log{hYo) after applying 
a mean internal extinction correction of 0.82 magnitudes. The gas cycling time of this sample is 
found to be igas = 7.5 ^2^ Gyr, and the volume-averaged equivalent width of the SINGG galaxies is 
EW{lla) = 28.8 ^47 A (21.2 Ig'g A without internal dust correction). As with similar surveys, these 
results imply that psfr{z) decreases drastically from z ~ 1.5 to the present. A comparison of the 
dynamical masses of the SINGG galaxies evaluated at their optical limits with their stellar and Hi 
masses shows significant evidence of downsizing: the most massive galaxies have a larger fraction of 
their mass locked up in stars compared with Hi, while the opposite is true for less massive galaxies. 
We show that the application of the Kennicutt star formation law to a galaxy having the median 
orbital time at the optical limit of this sample results in a star formation rate decay with cosmic time 
similar to that given by the psfr{z) evolution. This implies that the psFniz) evolution is primarily 
due to the secular evolution of galaxies, rather than interactions or mergers. This is consistent with 
the morphologies predominantly seen in the SINGG sample. 
Subject headings: galaxies: evolution - galaxies: ISM - galaxies: starburst - stars: formation - surveys 



1. INTRODUCTION 

The star formation rate density of the universe has 
changed considerably since 2; ~ 2, decreasing by ap- 
proximately an order of magnitu de. This decrease 
has been widely discussed (e.g. Madau et al. 1996; 
IPei. Fall, fc H'^^^s^^liggg: Somerville e.t nl. 2001 : Honkins 
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I2004J) because the evolution of the star formation rate 
density acts to constrain all models of galaxy forma- 
tion and evolution. Redshift-dependent luminosity den- 
sities (such as the R band luminosity density Ir{z)) 
and star formation rate densities {psfr{z)) remain some 
of the best constraints on these models. The value of 
Ir(z) at z Ki constrains the evolution of stellar mass 
IjMadau et oiJ 11998') . Likewise, the value of psfr(-z) at 
z sa helps determine the rel ative contributions of b urst 
and quiescent star formation l|Somerville et a0l200lD and 
the chemical evolution of the universe l)Pei fc Fal]l ll995'l. 
Estimates of psfr from previous studies span a factor 
of two or more, and do not always agree within their 
stated uncertainties, as shown in Table ^ AH surveys 
suffer biases, however, and these may explain the large 
discrepancies between the densities derived from each. 

The process of star formation leaves measurable sig- 
natures across the electromagnetic spectrum, allowing 
numerous methods for selecting star-forming galaxies, 
each with its own set of biases. For example, ob- 
jective prism surveys for emissio n line galaxies (e.g. 
iGallego et oiJll995llGronwalBl997!) result in a large, con- 
sistent bias towards galaxies with high surface bright- 
ness, high equivalent width emissio n lines (e.g. [Salzerl 



I1989D . Ultraviolet (UV) surveys (e.g. lTVever et alW U 
are biased against galaxies highly a ttenuated by dust. 
Conversely, far-infrared (FIR) (e.g. iFlo res et all 119991 
[Perez-Gonzalcz et al. 2005) or sub-millimeter surveys se- 
lect galaxies by the presence of dust-reprocessed opti- 
cal and UV light, resulting in a bias against galaxies 
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poss essing little dust. Si milarly, 1.4 GHz radio surveys 
(e.g. iSerieant et QiJl2002^ bias against systems with low 
dust content l|Bellll2003|) : in addition, the correlation 
between radio flux and FIR flux (which is then used 
to calculate star formation rate) is not linear, result- 
ing in a pronounced bias against low luminosity gal axies 
IjDeyereux fc EalesllT989t lYun. Reddy. fc Condonl l2001 ) . 

The 21-cm spectral line of neutral hydrogen (Hi) 
presents a useful starting point for star formation sur- 
veys. New stars form out of the interstellar medium, of 
which Hi is a key component. While it is the molecu- 
lar component of the ISM from which new stars form, 
a strong correlation has been found b etween Hi surfac e 
density and star formation intensity l|Kennicuttlll998j) . 
Hi flux is unaffected by dust extinction, and has been 
measured in every type of star-forming galaxy. As a re- 
sult, selecting targets by Hi mass creates a sample set 
free of the common optical biases. For this, we use the 
Hi Parkes All Sky Survey (HiPASS), an Hi survey of the 
southern sky [6 < +2°) over a velocity range from —1280 
to +127 00 km s^^ obta ined w ith the 21-cm multibeam 
receiver (iStave lcv-Smit h et a?jll9 96) at the 64m Parkes 
radio telescope (,Barnes et gi.ll200 1'l. 

We have commenced the Survey of Ionization in Neu- 
tral Gas Galaxies (SINGG) to survey star formation in 
the HiPASS galaxies. Meurer et al. (2006, hereafter Pa- 
per I) introduce the SINGG survey, including sample se- 
lection, methods used, and basic measurements. In total 
the SINGG sample set includes 468 Hi-selected targets 
covering the mass range 7.0 < log(A^Hi/A40) < 11.0 in 
approximately equal numbers per decade of mass. Dis- 
tance was also used in the selection process; the nearest 
sources were preferentially selected at each mass to im- 
prove physical resolution and avoid confusion. 

To measure star formation, SINGG uses Ha, the most 
readily accessible optical tracer of star formation at low 
redshifts. H« is a recombination line primarily result- 
ing from hydrogen photoionization. Because this re- 
quires ionizing UV photons, the majority of Ha flux will 
be produced near the most massive O stars, whose ex- 
tremely short lifespans (^10 Myr) make them good in- 
dicators of the current star formation rate. SINGG uses 
R band measurements for continuum subtraction, which 
also provide useful measurements of the existing stellar 
populations. Observations presented in Paper I consist of 
the 93 HiPASS extragalactic Hi targets fully processed 
to date. Due to the large beam size of HiPASS, 13 of 
these Hi targets contain between two and four distinct 
Ha objects. As a result, the SINGG sample includes a 
total of 111 individual Ha-emitting galaxies, which we 
refer to as SINGG Release 1, or SRI. We however ex- 
clude one Hi target, J0403-01, and its single Ha-emitting 
galaxy due to excessive sky uncertainties and foreground 
field contamination. The results included in this paper 
are entirely generated from the data presented in Paper 
I. 

This paper uses the SINGG observations to derive the 
Ha and R band luminosity densities of the local uni- 
verse, designated /hq(z) and Ir{z) respectively, where z 
indicates the mean redshift of the survey. For this and 
most of the other local (z « 0) studies, we will omit the 
(z) notation except where required. Combined with the 
Hi data from HiPASS, these yield the star formation 
rate density of the local universe, the stellar luminosity 



density, the density of the neutral ISM, and the cosmic 
gas cycling time. Section 12 explains the methodology 
used to determine the various volume densities and their 
uncertainties. Section |31 gives the results of our calcula- 
tions. Section 21 gives the results and compares to those 
of other surveys, while Section discusses some of the 
implications of our results. We use a ACDM cosmology 
{Qq = 0.3, r^A= 0.7) with a Hub ble constant of Hq — 70 
km s^i Mpc^\ and a lSalpeteJ {l955) IMF between 0.1 
and 100 Mq. 

2. METHODOLOGY 

2.1. Hi mass 

Since the SINGG sample is not volume complete, we tie 
our results to an Hi Mass Function (HiMF). Hi parame- 
ters are primarily derived from the data given in our two 
mai n source catalogs: HiCAT, the final HiPASS cata- 
log l|Mever fl< oiJl2^ . and BGC the H iPASS Bright 
Galaxy Catalog iJKoribalski et giJl2004D . These cata- 
logs have 95% completeness limits in velocity-integrated 
21 cm flux densities oi J Si, dv ^ 5 and 25 Jy km 
s^^, respectively. The Hi mass of a galaxy at a dis- 
tance of D Mpc is derived using the standard relation 
A^Hi [Mq] = 2.36 X lO^D^ j g^ j^ jjy ^^^ g-ij jRobertJ 
1962). As a result, the flux density limits for these two 
samples correspond to Hi mass limits of 1.2 and 5.9 
X 10^ D'^ M.Q, respectively. The SRI sources have dis- 
tances of 4 - 73 Mpc (with the majority falling within the 
10 - 20 Mpc range) and Hi masses of 10^-^ - 10^^-^ Mq. 
Detailed information about the Hi masses of the SINGG 
galaxies can be found in the original catalog publications 
and Paper I. 

Distances for most sources arc derived from radial ve- 
locities using the model of Mould et al. (2000), which 
corrects for infalls towards the Virgo cluster. Great At- 
tractor, and Shapley supercl uster. Distances to th e near- 
est galaxies are taken from iKarachentsev et al\ ((2004). 
This alters the distances and Hi masses for many 
SING G galaxies when c ompared to other HiPASS pa- 
pers IjZwaan et a^J 120051) . and marginally alters the Hi 
Mass Function used to correct for incompleteness. The 
values given throughout this paper are those derived us- 
ing the Mould model for the full HiPASS sample. 

We use the standard Schechter function to parameter- 
ize the HiMF: 
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where ^(A^hi) represents the number density of galaxies 
as a function of Hi mass (in Mpc^'^), Mi, is the char- 
acteristic mass, a is the "faint" end slope and 0* is the 
normalization factor. We c alculate Sche c hter fi ts to raw, 
binned 6{Miii) data from Zwaa n et al\ ll2005tl . and ad- 
justed to the distance of Mould ed al\ ()2000() for Hq = 
70km s^i Mpc-i. This yields M^, = 109-92±004 f^-2 
Mq, a = -1.41 ± 0.05 and 6^ = (3.86 ± 0.7) x lO'^ /ifg 
Mpc~'^ dex^^; these values are compared to other Hi 
Mass Functions in Table El As is usually the case, the 
errors on Mi, and 0* are highly correlated. 

For the 13 SRI Hi targets containing multiple Ho 
sources, HiPASS can only provide the total Hi mass 
for each target, with no ability to distinguish the con- 
tributions of each individual galaxy. When calculating 
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our luminosity densities, the luminosities of the individ- 
ual galaxies within each Hi target are combined and the 
total is treated as a single aggregate object. As the mass 
function given above was generated using similarly com- 
bined Hi masses, this approach should not substantially 
bias our results. 

2.2. Volume densities 

Using the HiMF, we derive the R continuum luminos- 
ity density Iji, Ho luminosity density /hqi and Hi mass 
density pni for the local universe. We denote values un- 
corrected for internal extinction with a prime (') symbol. 

The volume density of a quantity x is found using: 
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where x — Lua when calculating Zhq, x = Lji when cal- 
culating Ifi, X = 1.0 when calculating the number density 
n, or a; = TWhi when calculating pni- 

By separating our data into Hi mass bins and combin- 
ing the results within each bin, Eq. |21is replaced by: 
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Here, i represents the array of mass bins, while j rep- 
resents the individual HiPASS targets within each bin. 
Alog(A^Hi)i is the logarithmic width of each mass bin; 
our bins are 0.5 decades wide for all bins except the low- 
est (7.0 < \o%{MvaI Mq) < 8.0), as shown in Table|3 

The SRI Ho and R luminosities are given in Paper I, 
along with lists of targets containing multiple emission 
line galaxies. These luminosities are derived from fluxes 
extracted using elliptical apertures, supplemented with 
the flux from outer disk Hii regions. Full discussions of 
our R and H« flux extraction procedures are given in 
Paper I. 

To calculate the Star Formation Rate (SFR) for each 
galaxy, we ad opt the conver sion from Ha luminosity to 
SFR given bv lKennicutt et al. (1994). This relationship 
is derived using a singlepower l aw Initial Mass Function 
(IMF) having a Salpeti3 l|1955D slope and spanning the 
mass range of 0.1 to 100 Mq- The resulting conversion 

SFR [Me yr-^] = f ^H.[ergs^i] \ 

^ " -^ ^ V 1-26 X 1041 J ^ ' 

We adopt this conversion to maintain consistency with 
numerous other studies. The choice in IMF is a major 
source of systematic uncertainty, since the conversion of 
Lua to SFR depends entirely on the fraction of stars 
falling within the mass range responsible for ionizing hy- 
drogen. For example, the various IMFs of Scalo (198(1 1 
would decrease the con yersion factor to anyw here from 
1. 7x10^° to8.4x 10"° (|Kennicutt et g^JlTOQ^ . and that 
oflKroupal (1200 ID would give a conversion of 1.9 x lO''^ 
((Brinchmann et aL""2004). However, all Ha surveys suf- 
fer from this bias equally, and so as with Hq, we convert 
other Ha surveys to our chosen IMF before comparison. 

2.3. Flux corrections 

To calculate star formation rates, we first have to cor- 
rect the flux data for foreground and internal extinc- 
tion, [Nil] contamination, and stellar absorption. The 



foregro und extinction corrections are accomplished us- 
ing the iSchlegel. Finkbeiner fc Davia l|1998D extinction 
maps. Since our selection avoids the Galactic plane and 
the Magellanic clouds, the foreground extinction correc- 
tion increases the derived values of Zua and If/ by an 
average of only 10%. Stellar absorption, as explained in 
Paper I, causes our measurem ents to consistently under - 
estimate Ha flux by 2 - 6% (jBrinchmann et aZ.ll2004D . 
To adjust for this, we increase each of our Ha fluxes by 
4%. 

For the [Nil] correction and internal dust absorption, 
A(Ho)int, we adopt the relationships with R band ab- 
solute magnitude prior to dust absorption corrections, 
M'j^, given by Hclmboldt et al. (2004) and converted to 
the AB magnitude system, as explained in Paper I. Each 
galaxy's individual [Nil] /Ha correction is based on its 
Hi velocity, velocity width, and the narrow-band trans- 
mission profile, as detailed in Paper I. Integrated over 
our entire sample, the [Nil] correction decreases lua by 
15%; Ir is not affected. The internal extinction correc- 
tion A(Ha)int increases our estimate of Zhq by a fac- 
tor of 2.1 (0.82 mag). As noted in Paper I, we as- 
sume that the R band internal dust absorption, A{R)int, 
is half that of A(Ha)int, due to the well-known phe- 
nomenon of increased extinction in Hii regions com- 
pared to the field (Fan clli. O'Conneh. & Thuan 198^ 
iGalzetti. Kinnev fc Stor chi-Bcrgman 1994), and so Ir in- 
creases by 0.41 magnitudes after extinction correction. 

3. RESULTS 

3.1. Volume densities 

The R band luminosity density, corrected only for 
Galactic extinction, is found to be 

I'r = (4.4 ± 0.7) x 10^^ hjQ erg s^^ A"^ Mpc'^ (5) 

while the Ha luminosity density, corrected for [Nil] con- 
tamination and Galactic extinction, is 
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(9.4 ± 1.8) X 10^^ hro erg s^^ Mpc" 



(6) 



With our adopted IMF and resultant star formation rate 
conversion, Eq. ^ the local star formation rate density 
(uncorrected for internal extinction) is 

fog(PsFR [-^0 yr~' Mpc-3]) = -2.13 t"o'oliran.) (7) 

±0.03(sj/s.) + log(/i7o) 

Corrected for internal dust extinction, these equations 
become: 

Iji =7.0 to,l X 10^^ /i7o erg s^^ A~^ Mpc'^ (8) 

fee - 2.0 toi X 10^9 hjo erg s-^ Mpc-^ (9) 

fog(psFR [Mq yr-iMpc-3]) = _i.80 ig;i?(mn.) (10) 

±Om{sys.) + log{hn) 

The uncertainties in each variable are explained in detail 
in subsections 13.21 and 13.31 
Additionally, we derive the Hi mass density to be 

PHI = 5.17 ± 0.38 X 10^ hro Mq Mpc"^ (11) 

and the mean number density of Hi-rich galaxies to be 

n - 0.112 t^oll hro Mpc'^ (12) 
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within our mass range, 7.0 < \og{MBil Mq) < H-O- 
Unlike the luminosity densities, this number density is 
completely dependent on the lower boundary chosen, 
as the integrated function increases as A^hi decreases. 
Fig- n shows the dependence of the luminosity densities 
(and, by extension, derived quantities such as psfr) on 
A^Hi, along with the systematic uncertainty in the to- 
tal value due to each bin. The plotted quantity is the 
fraction of the total luminosity density coming from one 
decade of mass. 

For comparison, we also derive psFR and pni for other 
mass functions, using the simpler distance model favored 
by HiPASS; these values are given in Table |21 Only 
the HiMF weighting of each Ho source is altered; the 
distances and fluxes remain unchanged. While the dif- 
ferences in individual HiMF parameters are relatively 
small, the resulting integrated mass and luminosity den- 
sities in Table|21vary by 0.20 - 0.22 dex, primarily due to 
the low-mass slope, a. While this discrepancy is compa- 
rable to the total uncertainties in each density, it is signif- 
icantly larger than the ^ 0.03 dex uncertainty caused by 
the uncertainties in the HiMF parameters themselves. 

3.2. Random errors 

Since the calculations of luminosity densities involve 
many different variables, most of which have their own 
uncertainties, the simple methods of error propagation 
would result in heavily correlated uncertainties. To bet- 
ter quantify the random uncertainties, we utilize Monte 
Carlo and "bootstrap" algorithms. 

We consider the following sources of random error: (1) 
the HiMF: the limited number of HiPASS sources used 
to derive each HiMF results in uncertainties in the pa- 
rameters used to fit that HiMF, and hence on luminos- 
ity densities; (2) error due to the limited SINGG SRI 
sample; (3) R and Ha flux uncertainties due to sky sub- 
traction; (4) Ho flux uncertainty due to continuum sub- 
traction; (5) the uncertainty in the Ho flux calibration, 
and (6-7) error due to the dispersion in the fits used to 
generate our [Nil] and ^(Ho)int corrections. The random 
errors are presented in Tabled we consider each of these 
terms in detail. 

To estimate the uncertainty due to the HiMF, we cre- 
ate one hundred realizations of the HiMF using a boot- 
strap resampling of the original data in HiCAT and 
adopting the Mould et al. (2000) distance model. The 
HiMF in each realization is created using the same two 
dimensional stepwise maximum likelihood technique of 
iZwaan et al\ (|2005). Each reahzed HiMF is fit to a 
Schechter function, which is used in Eq. |31 to determine 
the resulting luminosity densities of each realization. The 
random uncertainty in each luminosity density due to the 
HiMF is then the dispersion about the mean luminosity 
density for all the realizations. 

To quantify the SRI sampling error we use a "boot- 
strap" resampling method, drawing 92 objects at ran- 
dom (with duplication allowed) from our data set. As 
with our other uncertainties, this randomization is re- 
peated ten thousand times, with the overall sampling 
uncertainty defined as the standard deviation of the re- 
sulting distribution. 

The errors due to the uncertainties in the continuum 
scaling ratio and background sky level for each galaxy 
are quantified through another form of Monte Carlo logic. 



For each of ten thousand iterations, we vary the sky level 
or continuum ratio within Gaussian distributions having 
the uncertainties derived from our error models, which 
in turn alter the measured flux (and by extension, the 
star formation rate) for each galaxy. Again, the lumi- 
nosity densities are recalculated for each iteration, and 
the dispersion about the mean is quoted as the resulting 
uncertainty. Additionally, we have an uncertainty due to 
our flux calibration method, as explained in Paper I; we 
estimate this uncertainty to be 0.04 magnitudes for im- 
ages requiring our 6568/28 narrow-band filter, and 0.02 
magnitudes for all others. 

To find the random error due to our A(Ho)int and [Nil] 
correc tions, we use the M'^ fits given in lHelmboldt et all 
l)2004() . each of which has a dispersion of 0.23 dex (Helm- 
boldt, priv. comm.). We make a series of realizations 
of our sample; within each realization, each galaxy's 
log(A(Ho)int) or log(i^[Nii]6583/^Ha) Correction is per- 
turbed by a Gaussian random deviate with the above 
dispersion. Luminosity densities are re-derived, and the 
dispersion about the mean simulated luminosity density 
is our random error. Separate sets of realizations are 
done to determine the errors due to A(Ho)int and those 
due to the [Nil] correction, with all other terms held 
fixed. 

Applied to the SINGG sample, we find the uncertain- 
ties as listed in Table 01 We find that the random uncer- 
tainties are dominated by the internal dust extinction 
(for the corrected luminosity densities only), the SRI 
sampling error, and to a lesser extent the HiMF uncer- 
tainties. While the sky and continuum subtraction un- 
certainties dominate the measurements of many individ- 
ual galaxies, when evaluated over our entire sample their 
contributions to the error budget are relatively small. 

3.3. Systematic errors 

The SINGG results also suffer from a series of system- 
atic uncertainties. We consider the following sources of 
systematic error: (1 - 2) uncertainty in the zeropoints 
of our [Nil] and 74(Ho)int corrections; (3) variation due 
to our choice of distance model; (4) variation due to our 
choice of HiMF. 

In addition to the dispersion mentioned in Section [3. 21 
our A(Ho)int a nd [Nil] corrections inclu de an uncertainty 
in the fit from iHelmboldt et al\ l|2004D itself. This zero- 
point uncertainty corresponds to the discrepancy in lu- 
minosity density for fits one standard deviation of mean 
away from the best fit, generated from the 196 sources 
of Janscn (2000). 

Our distance model is a significant source of systematic 
uncertainty, as most SINGG targets are located within a 
distance of 20 Mpc. Variations in the distances used can 
result in large changes in observed luminosities as well as 
the underlying Hi Mass Function, and as a result alter 
the derived star formation rates. To quantify this effect, 
we ca lculate psfr using both our default ^ould et al\ 
(|2000) distanc e model as well as t he simpler local-group 
model used bv IZwaan et a,l\ l|200fi() . and quote the differ- 
ence as our uncertainty. 

Finally, we calculate our luminosity densities using a 
variety of Hi mass functions and the local-group distance 
model, with results shown in TablcEJ While the HiPASS 
team has produced several Hi mass functi ons, each su- 
persedes the one before, with the work of IZwaan et al\ 
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l|2005(l comprising the most complete version. As a re- 
sult, while the resulting values of psFR vary by non- 
negligible amounts, the earlier HiMF values are only for 
comparison purposes, and we do not include this error in 
our final uncertainties. 

Our quoted systematic uncertainties are completely 
dominated by our choice of distance model. 

3.4. Equivalent width 

We define the volume-averaged Ha equivalent width to 
be the ratio of Ha flux density to R band flux density. We 
flnd the volume-averaged equivalent width of our sample 
to be 

EW'iRa) = 21.2 tti A (13) 

without internal extinction corrections, and 

EWiKa) = 28.8 tlj A (14) 

after all corrections have been applied. The difference 
between these two values is due to the differential na- 
ture (line versus contin uum) of the extinct ion law applied 
IjCalzetti. Kinnev fc Storchi-Bergman 1994, Paper I). 

The Ha equivalent width is a measure of the rate of 
transformation of the ISM into stars, relative to the 
existing stellar population. Comparison of the dust- 
corrected equivalent w idth to the theoretical curves of 
iKennic^itt et al\ (|1994fl yields the birthrate parameter b, 
corresponding to the ratio of the star formation rate at 
the present time to that averaged over the age of the 
galactic disk. This interpretation is heavil y dependent 
on the IMF used; for our adopted ISabeteil l[l955. ') IMF, 
we find 



b = 0.26 



ho.io 



(15) 



which supports the common assertion that the cosmic 
star formation density has decreased substantially since 
z ~ 1.5 within star- forming galaxies. Other IMFs result 
in greater values of 6, ranging as high as 2.8 for bottom- 
heavy IMFs such as lScaloll|1986(l . 



3.5. Gas cycling timescale 

Using the SINGG data, we determine the gas cycling 
timescale igas, the time it would take for star formation 
to process the existing neutral and molecular ISM of a 
galaxy at its current rate of star formation. This value 
can be derived from Hi data if two assumptions are made. 
First, we assume that Hi emission occurs predominantly 
at low optical depth. Second, we assume that the ratio of 
molecular hydrogen to neutral atomic hydrogen remains 
constant, independent of the other characteristics of the 
galaxy. 

To calculate igas we use the equation given in Paper I, 
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where Ma 



2.3 A^HT, derived from the observations 



of lYoung et al\ ()1996|) , accounts for the typical molecular 
hydrogen as well as the helium content of the ISM. 

The values of tgas for our array of mass bins are given 
in Table 13 The volume-averaged gas cycling timescale 
for our sample can be derived from tgas ~ 2.3 (phi/psfr); 
with the values for pni smd pspR presented in Section f3.1l 
we find this timescale to be igas — 7.5 Ij ^ ^y^~- This is 
somewhat less than the Hubble time, consistent with the 



previous findings of lKennicutt et all ^1994) . Fig.|3shows 
the gas cycling timescales for each of the SINGG galaxies, 
as well as the average within each Hi mass bin. While 
the galaxy-to-galaxy variation in igas is substantial (0.54 
dex), the mean value does not appear to vary greatly 
with A^Hi- 

Our choices of distance model and Hi Mass Function 
have much smaller impacts on the cosmic igas than they 
do on psFR) and have no effect on the gas cycling times 
of the individual galaxies within our sample. Using the 
other Hi Mass Functions given in Table[21result in cosmic 
tgas values of 7.0 to 7.3 Gyr, well within the error bars 
of our adopted HiMF. 

One should not take these gas cycling timescales too 
literally, as several of the assumptions behind them may 
be questioned. Much of the neutral interstellar medium 
may have a large optical depth, which would imply larger 
values of A^hi than measured, causing fgas to be some- 
what underestimated. This effect probably results in 
a dis crepancy of less than 20% IjHavnes fc GiovanelH 
Il984( ). As mentioned in Paper I, the observed ratio of 
(A^H2/A4hi) used to derive Eq. E|has a dispersion of 
0.58 dex, or almost a factor of 4; the ratio of CO lumi- 
nosity to H2 mass also has a large uncertainty. Hence, 
the tgas estimate of any single galaxy is likely to be highly 
uncertain. Similarly, we have not accounted for system- 
atic effects. It is well known that more massive and 
higher surface brightness galaxies are easier to detect in 
CO emission than dwarfs. In addition, our scenario does 
not include the "hot phase" of the ISM. This accounts 
for the X-ray emitting halos around galaxies as well as 
the intracluster medium and the intergalactic medium. 
This is li kely to be the largest bary onic component of the 
universe IjFukugita. Hogan. fc Pee bles 1993- We also do 
not account for the return to the ISM of material pro- 
cessed by stars, the stellar yield. Because of these limita- 
tions we are careful not to imply that tgas is a consump- 
tion timescale; rather, it is simply the time it would take 
for the present rate of star formation to form the mass 
of observed neutral ISM and inferred molecular ISM into 
new stars. 

3.6. Dynamical parameters 

The dynamical mass located within the optical radius 
of a galaxy, A^dyn, is a useful way to estimate its mass 
(including both baryonic and dark matter) from easily 
observable quantities. The corresponding orbital time, 
torb, is also used for comparison with various models 
of galaxy evolution. To flnd these values, we flrst ap- 
proximate the circular velocity of the edge of the galaxy 
from the FWHM spread in the velocities measured by 
HiPASS, corrected for the estimated inclination of each 
galaxy. That is. 



W; 



50 



^2 sin{i) _ 
where our inclination angle, i, is deflned as 



sin (i) 



(17) 



(18) 



lIBottinelli. Gouguenheim. Patural. fc de VaucouleursI 
Il983|) . where a.c are the axial lengths of the oblate 
spheroid fit to each galaxy, and a, b are the axial lengths 
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of the elliptical aperture used to measure optical fluxes. 
Since we do not wish to include any dependence on mor- 
phological classifications in the SINGG methodology, we 
assume a constant {c/a) = 0.20, except in three cases 
where our estimated {h/a) is less than this amount, 
where we simply set c — b. To reduce uncertainty 
in TWdyn we exclude the 13 Hi targets with multiple 
emission line galaxies, as well as 14 additional face-on 
galaxies {{b/a) > 0.8, or i < 38°). 

We then find the mass contained within the observed 
limit of each galaxy to be 



A^dyn = 



and the orbital time to be 



^circ ' max 



G 



^orb — 



27r Tn 



^cir 



(19) 



(20) 



A^dyn is evaluated using the radius at which we trun- 
cate our light profiles, Tmax- This radius is found using 
only the optical light of the galaxy, therefore matter lo- 
cated outside it will not be included. A^dyn is not a true 
total mass; rather, it is a crude estimate of the mass 
contained within the optical radius. Since Hi typically 



extends beyond 



with a flat rotation curve, the total 



mass (including dark matter) is larger, and our A^dyn 
estimates are thus lower limits to the total mass. A^dyn 
is compared with our HiPASS-derived values of A^hi in 
Fig. 13^, and with our observed R band luminosity, L^, 
in Fig. 1213. After excluding galaxies as explained above 
and fitting with a robust least absolute deviation fit, we 
find: 

log(A^dyn) = 1.26 log(7WHi) - 1-36 ± 0.37 (21) 
log(7Udyn) = 0.79 log(Li^) + 2.83 ± 0.36 (22) 

Eq. [21 implies that the fraction of A^dyn made of Hi is 
typically more than six times higher for our lowest Hi 
masses (A^hi ~ 10^'^ -^e) than it is for those at our 
high-mass extreme (A^hi ~ lO"'^'''^ -Mq)- Likewise, if 
we were to assume that A4/L remains constant over our 
range of dynamical masses. Eg. 1221 would imply that the 
visible mass fraction is up to 10 times larger for our most 
massive galaxies as for our least massive. For any fur- 
ther equations requiring A^dyn: we use the fit of Eq. TF3 
to set the dynamical mass of those galaxies exluded from 
our fit because of inclination or multiple Ha sources. 
While Eq. [221 is reminiscent of the classic relations of 
iTullv fc Fisher l)1977j) . our A^dyn depends both on the 
width of the Hi velocity profile and the optical radius, 
'"max- As a result, our result should not be directly com- 
pared to the TuUy-Fisher relationship, Mi, ~ ^circ- 

We are primarily interested in torb because it lets us 
address the redshift evolu tion of SFR in con junction with 
the star formation law of'Kennic uttI l|1998l) . which states 
that the global star formation law of galaxies is equiva- 
lent to galaxies converting 21% of their ISM mass within 
one orbital time iorb evaluated at the radius where the 
Hii region distribution is truncated, rnii- We have not 
measured run as such, for our sample. However, Paper 
I provides two radii that should bracket this: r9o(Ha), 
the radius enclosing 90% of the Ho flux, and r„iax; for 
most sources, rmax contains all the discernable emission 
in both R and Ha. We evaluate the orbital time, torb, 



at these two radii. As with our determination of A^dyn, 
iorb is ill-determined for galaxies which have i < 38° or 
are in multiple ELG systems. In those cases we estimate 
log(torb) from Lji using least absolute deviation fits, sim- 
ilar to Eq. E21 

log(iorb(rmax)) =0.099 log(Lfl) - 1.04 ± 0.20(23) 
log(iorb(r9o(Ha))) -0.106 log(Lfl) - 1.35 ± 0.23(24) 

These fits are then used to set orbital times for the ex- 
cluded galaxies, in the same manner as for Aldyn- 

As with Hi mass, we use TWdyn to find the local dy- 
namical mass density, Pdyn, with Eq. O We find 



Pdyn 



9.3 tli X 10« 



/l70 Mq Mpc 



(25) 



Again, we emphasize that this is only the mass density 
of the local universe residing within the optical radii of 
Hi-rich galaxies. Of this density, 15% is found in galax- 
ies with multiple Ha sources, while an additional 34% 
is found in the face-on galaxies excluded from Eq. 1221 
Comparing to our Hi mass density, pni, we estimate that 
5.6% of the dynamical mass of local galaxies consists of 
neutral hydrogen in stars. 

4. DISCUSSION 
4.1. Completeness 

One of our goals is to determine how representative an 
Hi-sclccted sample is. To do this, we compare the R band 
luminosity density Z^ and dynamical mass density pdyn 
derived from an Hi-selected sample to those derived from 
more "complete" samples, such as the full SDSS sample 
set. 

Our uncorrected value for I'j^, 4.4 ± 0.7 x lO'^^ h^o erg 
s~^ A^^ Mpc~'^, compa res well to the SPSS-d erived ^-^r 
and '^■^i band values of iBlanton et al\ l)2fl03D . 6.17 and 
6.70 X 10'^^ hiQ erg s~-^ A~-^ Mpc~^ respectively. Inter- 
polating between these two bands by wavelength gives 
I'ji « 6.30 X 10^7 hrp erg s"! A"! Mpc'^ for the Sloan 
survey, which we will use as our "cosmic" R band den- 
sity. Our estimated Vj^ is approximately 70% of this 
value; therefore, the majority of the stars comprising 
the derived SDSS R band luminosity density are located 
in galaxies containing measurable quantities of Hi. If 
we were to assume that our I'j^ is low due to gas-poor 
galaxies (ellipticals) being absent from our sample, while 
our I'yi^ is complete, then we can derive a corrected 
"cosmic" EW{ila) « lHa/{lR/0.70) = 20.3 tH A. 
This reduces the Salpeter-derived birthrate parameter to 
u ~ u.io _o.o3- 

Note that the SRI sample size implies that certain 
rare types of galaxies simply won't be represented in 
these results. For instance, ultral uminous IR galaxies 
(ULIRGs, defined as Lir > 1 0^^ T.^. ls^ers fc Mirah^ 
il596) are rare enough that the volume contained within 
the HiPASS redshift and declination limits would only 
include roughly one ULIRG. As the sample used in this 
paper includes only 92 out of the 4315 HiPASS targets, 
it was extremely unlikely that any ULIRGs would be in- 
cluded in our observations. 

When our dynamical mass density pdyn is expressed 
as a fraction of the critical Einstein-de Sitter mass den- 
sity {pcrtt = 1-36 X 10" /i^o Mq Mpc~3), we find 
^dyn = 0.0068 hjQ, or 2.3% of Qq in the concordance 
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cosmology. For comparison, iCole et M l)200H) estimates 
the the z « mass density of stars, ristars, to be between 
0.0029 (from an IR luminosity function) and 0.0020 (from 
st ar formation tracers ). Similarly, the Hi Mass Function 
of lZwaan et al\ (|200,'? ) gives the local Hi mass density to 
be « 0.0004 hjQ, which would imply a gas density on the 
orderof rigas = 0.0009 h^Q using the same logic as in Sec- 
tion |^3 As a result, our pdyn implies that around half 
of the mass located within the optical disks of nearby 
galaxies consists of gas and stars, with the remainder 
most likely consisting of dark matter. 

This comparison must be viewed with some caution; 
our pdyn is a crude estimate of all matter within the op- 
tical radii of Hi-selected galaxies, including substantial 
quantities of dark matter. However, it excludes galac- 
tic gas (especially Hi) extending beyond the optically- 
selected Tniax, and also samples little of the hot plasma 
which resides in galaxy clusters and the intergalactic 
medium. We also do not sample any of the mass in Hi- 
free galaxies, predominantly early-type galaxies. This 
comparison only shows that the mass we sample is com- 
parable to the baryon content of the local star-forming 
galaxies, and that the ratio of baryonic matter to dark 
matter within these galaxies is substantially different 
than that of the universe as a whole. 

4.2. psFR and lua in context 

To compare our value of p s fr to those found by other 
studies, we refer to IHopkina IJ2004J) . which compiled the 
results of 33 other star formation rate density p apers and 
correc ted each to a uniform ACDM cosmology, a lSalpeten 
l|1955D IMF, and with a Hubble constant of Ho = 70 km 
s~^ Mpc"\ f7o == 0.3, and Q.k= 0.7. We also include th e 



IR-derived psfr(-z) data ofl Perez- Gonzalez et al\ ()200£ 
The results before internal extinction correction are plot- 
ted in Fig. 2ti while the corrected values are plotted in 
Fig.ElD. 

In Paper I we assert that SINGG is inherently less 
biased than other surveys, due to selecting sources by 
ISM content at radio wavelengths. Objective-prism sur- 
veys limit their samples by equivalent width and surface 
brightness, while an Hi-selected sample can include dif- 
fuse sources or those with little star formation. If this 
is the case, we should expect to recover more star for- 
mation in the local universe. However, there are many 
steps involved in turning measured fluxes into lumi- 
nosity density estimates; since techniques and assump- 
tions vary between groups it is important to compare 
our results in as consistent a form as possible. Few 
previous studies present actual Ho luminosity densities, 
and so our primary comparison will be between values 
of /5gFR(^) corrected to our adopted IMF and extrapo- 
lated to z = 0. While some studies have only provided 
extinction-corrected psfk{z) estimates, in most cases we 
are able to work backwards using published corrections 
to derive /o'sfrI-^) foi' each. 

When the surveys located at 0.0 < z < 1.0 are linearly 
fit to a simple P's,fk[z) — psfr(O) (1 + zY relationship, 
we find (3 = 3 .00 ± 0.13, which mat ches well with the 
13 = 3. 2 til oflT^ y\nc^\\ et al\ IM)^ and the /3 = 3.1 ± 
0.5 of iPerez-Gonzalez et al\ (2005). As a result, we will 
use this proportionality to extrapolate the star formation 
rate densities of low-redshift surveys to z ~ 0; the results 



are given in Tabled 

The relative positions of the data points in Fig. ^,b 
change due to the internal extinction correction. While 
the spread in points noticeably decreases between the 
two figures, we cannot claim to measure the highest cor- 
rected psfr(O). This is primarily due to our relatively 
mild extinction correction of 0.82 magnitudes, while the 
other samples noted in Table ^ have corrections rang- 
ing from 1.0 to 1.4 mag. Although mild, our correc- 
tion correlates well with unpublished work from the 2dF 
Galaxy Redshift Survey (Folkes et al. 1999) . which used 
the Balmer decrement to estimate the Ha extinction of 
160,000 line emitters; the typical value was found to 
be approximately 0.8 magnitudes. This discrepancy in 
^mt(Ha) is most likely caused by the same luminosity 
bias mentioned above. 

Even with a value of psFR comparable to other surveys, 
it is too early to declare a consensus because of the differ- 
ences between the samples. For instance, we assert that 
the objective-prism selection method results in a large, 
consistent bias against non-starburst galaxies. Here we 
define starburst galaxies as those with Ha equivalent 
widths within their half-light radii of EW^Q{Ra) > 50A, 
without c o rrecti on for internal extinction. According to 
iHeckmaa ()1998j) . starbursts are estimated to comprise 
15 - 20% of the population of the local universe. In the 
UCM spectroscopic survey, 72% of the sample consists 
of galaxies with equivalent widths above this threshold. 
For SINGG SRI, 16 out of the 110 Ha sources (14.5%) 
met this criterion; these sources are collectively respon- 
sible for 25% of our final star formation rate density. 
If we assume that the Hi-selected SINGG sample set is 
not significantly biased towards or against these galax- 
ies, emission line-selected surveys (such as UCM) should 
be significantly underestimating the value of psfr due to 
underrepresentation of the non-starburst galaxies which 
appear to produce the majority of the star formation in 
the local universe. 

4.3. Breakdown of luminosity density and psFR 

Fig. El shows the observed contributions to the R band 
and Ha luminosity densities as a function of a variety of 
different quantities; this allows a comparison between our 
sample and other local samples of galaxies. Table^lgives 
the values of each parameter at the 10*'*, 25*'', 50"*, 75*'', 
and 90*'' percentiles of each luminosity density. In all 
cases, extinction and [Nil] corrections have been applied 
as appropriate. 

Plot (a) relates the luminosity densities to the HiPASS 
Hi mass. We find that 72% of Ir and 70% of lua are 
found in galaxies with Hi masses below our value of 
7W*(Hi) . However, 77% of pni falls below this value, 
so the stellar luminosity density and psFR are slightly 
weighted towards higher Hi masses than pni- 

Plot (b) relates the luminosity densities to the dynam- 
ical mass, with the fit from Fig. ^p used to estimate dy- 
namical masses for those galaxies with multiple sources 
or low axial ratios. Galaxies with low dynamical masses 
contribute substantially more to the Ha luminosity than 
the R band density, and the SINGG sample extends 
across a wide range of A^dyn values. 

Plot (c) relates I to the R band luminosity by way of 
the absolute magnitude, M^. As expected, galaxies with 
high R band luminosities contribute more to Ijj than to 
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fart ■ For comparison, the study of iBrinchmann et al\ 
l^20^)4^ used a stellar mass function with its knee at 
Mi, = 10^^-^^ Mq] as this is a stellar mass (not simply Hi 
mass), it should be approximately proportional to the R 
band luminosity. The 10* ^ and 50*'* percentiles for fa^ in 
IBrinchmann et al\ ((20^ were quoted as log(7W^) — 9.0 
and 10.3, respectively, a difference of 1.3 dex. For com- 
parison, the corresponding SINGG values of M^j for those 
percentiles differ by 4.1 magnitudes (1.7 dex), as shown 
in Table 

Plot (d) relates I to the star formation rate of each 
galaxy (directly proportional to Ho luminosity, as shown 
in Eq. 01) As expected, galaxies with high Ha luminosi- 
ties contribute more to Zhq than Z^. However, the frac- 
tion of Ha luminosity density caused by luminous galax- 
ies in SINGG is significantly lower than in other surveys. 
Only 12.2% of the SINGG faa is caused by galaxies with 
observed star formation rates greater than 10.0 A4q yr~^ 
once extinction corrections have been applied. For com- 
parison, we integrate the luminosity functions quoted by 
other surveys; the UCM study of Galleao et al. ( 19951 
has 26.5 % of its fa„ come from gala xies above this thresh- 
old, while lTresse fc Maddoxllll99^ has 23.1%. At the ex- 
treme cases. iSullivan et a? j 1120001) has only 9.3%, while 
iPerez-Gonzalez et al\ 1)2003") has 51.4% of its luminos- 
ity density come from galaxies with star formation rates 
above 10.0 M.q yr~^. This further illustrates that the 
apparent consensus in psFR may hide important differ- 
ences in the various surveys of local star formation. 

Plot (e) relates I to the narrow-band internal dust ex- 
tinction correction, A(Ha). While the extinction correc- 
tion in our overall Ha luminosity density is 0.82 magni- 
tudes, only 37% of the R density comes from galaxies 
with smaller values, compared to 47% of the Ha den- 
sity. That is, Iji is weighted towards the more luminous 
galaxies, which have higher internal extinctions. 

Plots (f ) and (g) relate I to the effective (50% flux) radii 
found in the R and Ha images, while plot (h) relates / to 
the ratio of these two radii. Again, these plots support 
the "downsizing" model of star formation; galaxies with 
large radii tend to contribute more to the R band lumi- 
nosity density than the Ha density, while galaxies with 
small radii contribute substantially more to Zua than 
In- The luminosity density contributions are evenly dis- 
tributed over our range of (re(Ha)/re(i?)). Galaxies with 
centrally-located star formation (re(Ha)/re(i?) < 1.0, 
suggestive of starbursts) do not dominate the luminos- 
ity densities; this contradicts what has be en observed in 
other studies, such as the UCM survey of iGalleeo et al\ 
Ip95l) . 

In all of the above plots, the smaller, less massive galax- 
ies contribute substantially more to the SINGG Ha lu- 
minosity density (and by extension, psfr) than they do 
to the R band density. This reinforces the "downsizing" 
model, where star formation activity shifts to smaller 
galaxies over time (Cowic. Songaila. Hu. & Cohen 1996). 

Plot (i) relates I to the ratio of the 90% flux ra- 
dius to the half-light radius in the R band. The in- 
verse of this ratio (referred to as the "concen t ration in- 
dex") has been examined bv lShimasaku et al\ l)2001t ) us- 
ing SDSS data. In that work, it was estimated that 
the boundary between bulge-dominated "early-type" and 
disk-dominated "late-type" galaxies occurred at a ratio 



of rgo (r')/re(r') = 3.03. Given that definition, we esti- 
mate that 71% of our R density ((3.2 ± 0.5) x lO^^/iyo erg 
s"-'^ A~^ Mpc~^ before extinction correction) and 82% of 
our Ho density can be attributed to "late-type" galax- 
ies, assuming the r-derived boundary does not change 
when shifted to our R band. For comparison. iHogg et al\ 
((2OOI find that 38% of the SDSS ^-^i band luminosity 
density comes from "red-type" (assumed to be bulge- 
dominated early-type) galaxies, meaning that a total of 
4.1 X 10^'^ /i7o erg s~^ A~^ Mpc~^ can be attributed 
to late-type galaxies. As their definition of red galaxies 
was conservative, this percen tage should include some 
early-type galaxies. Likewise, iBal drv et al\ l)2004f l fit a 
bimodal distribution to observed galaxy colors to deter- 
mine that 58% of the SDSS "'^r band luminosity density 
comes from blue galaxies (for a density of 3.6 x 10^^ ft.70 
erg s"^ A^i Mpc"3). As a result, while the SINGG sur- 
vey recovers around 70% of the SDSS In, we recover a 
larger fraction (80 - 90%) of the luminosity density from 
late-type galaxies. 

Plot (j) shows the dependence of I on Ho equiva- 
lent width. T he /Jrfr derived from the UCM spectro- 
scopic survey ()Gallego et al\ \19i^^ uses only galaxies 
with equivalent widths larger than 10 A. According to 
our data, 24.1% of the local extinction-corrected R band 
luminosity density and 4.5% of the corrected Ho lumi- 
nosity density come from galaxies with equivalent widths 
below this threshold. 

Plots (k) and (1) relate I to surface brightnesses, fJ,e{R) 
and Se(Ho). As expected, galaxies with high Ho surface 
brightnesses (i.e., starbursts) contribute a much larger 
fraction of lua tha n oi Ir. F or comparison, the SDSS- 
derived sample of Blan ton et a,l\ (|200 5,1 has a surface 
brightness limit of /ie(_R) = 24.0 ABmag arcsec"^; 1.1% 
of our In and 1.2% of /hq come from galaxies below this 
cutoff. 

4.4. Recent redshift evolution in y6sFR(^) 

Fig. ^ shows that the psfr{z) has declined by ap- 
proximately a factor of ten from z ~ 1 to the present. 
Specifically, fitting the internal dust corrected data with 
z < 1 with a robust (outlier resistant) linear fit 



log(psFR(^)) = log(psFR) + vz 



(26) 



yields log(psFR [Mq yr"! Mpc'^]) = -I.73 and 77 = 
1.02 dex per redshift. We can estimate the systematic 
uncertainties in this fit by categorizing the psFniz) data 
by the star formation tracer used in each measurement: 
(1) optical emission lines; (2) ultraviolet continuum; and 
(3) FIR or radio continuum. This yields log(pgFR) — 
-1.80, -1.53, -1.94, and ry= 1.05, 0.43, 1.37 for survey 
types 1, 2, and 3, respectively. This suggests that the 
uncertainties in the zeropoint and slope are 0.20 and 0.47 
dex, respectively. 

What c auses the decrease i n Pfiv n(z) with cosmic time? 
Recently, iKauffmann et all l)2004f) presented a model of 
the recent star formation history of galaxies in the local 
universe that provides a useful context to answering this 
question. In their model, star formation events occur 
when galaxies merge; subsequently, the SFR inside the 
galaxy decays according to a pres cription sim ilar in na- 
ture to the star format ion laws of lKennicuttI 11998) and 
iDopita fc Rvdeill|2005j) . They use a high-resolution Cold 
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Dark Matter (CDM) simulation, identify galaxy merg- 
ers with CDM halo mergers, and parameterize the form 
of their star formation law by the average stellar mass 
density. Essentially, the physics of their model can be 
separated into CDM effects (the merger of halos), and 
baryonic physics (the star formation law). We can then 
rephrase the question above: is the decrease in psfr(^) 
with cosmic time driven by CDM or baryonic physics? 
Are we seeing a recent decrease in the merger rate of ha- 
los? Or, are we seeing the secular decay of the SFR after 
a much earlier decrease in the halo merger rate? 

The dynamical information on our sample provided by 
the Hi line widths provides a means to address this ques- 
tion. The orbital time, iorb, has been evaluated at two 
radii, rgo(Ha) and Tmax, as explained in Section [3. 61 Us- 
ing the same methods as those in Section lOl we find that 
the 25%, median, and 75% contributions to lua (and, by 
extension, psfr) in our sample occur for torb = 0.30, 0.57, 
0.60 Gyr when evaluated at rgo{iia) and 0.56, 1.00, 1.03 
when evaluated at rmax- 

We can use the orbital times to estimate the change in 
SFR with redshift for the galaxies in our sample obeying 
Eg. 1261 If we assume that the star formation law remains 
invariant over this timescale, then consumption of 21% 
of the ISM is equivalent to a 21% decrease in the SFR 
over iorb (Kcnnicutt 1998). Translating iorb to a look 
back time, the logarithmic change in SFR is given by 



<51og(SFR[Xo yr-i]) 1.43 /i 
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in the same units as the slope rj in Eq. 1261 The median 
iorb translates to (51og(SFR)/5z = 2.5, 1.4 as defined by 
r9o(Ho) and rmax, while the interquartile ranges are 4.8 
- 2.4 and 2.6 - 1.4, respectively. We see that the SFR of 
the galaxies in our sample will decay over a wide range 
of timescales. However, a galaxy with iorb similar to the 
median galaxy in our sample will have its star formation 
decay at a rate similar to that observed in the universe 
as a whole. 

5. CONCLUSIONS 

Our measurement of the local star formation rate den- 
sity, log(psFR [Mq yr-1 Mpc"^]) = -1.80 +log(/i7o), is 
similar to those of most previous z ^ studies, as shown 
in Table Q] When compared to the results of other sur- 



veys, this agrees with the consensus that the star forma- 
tion rate density has decreased by an order of magnitude 
from a redshift of z ~ 1 to the present. While a consen- 
sus in psFR appears to be emerging, this may be largely 
illusory. There remain significant differences between the 
various surveys in terms of extinction corrections and the 
contributions of various types of galaxies to the totals. 
Systematic effects are large and need to be accounted for 
before a true consensus can emerge. 

The decreasing psfr{z) with cosmic time has impor- 
tant implications for models of galaxy evolution, as this 
trend suggests a shift from the faster "burst" star for- 
mation process to the slower "quiescent" process, with 
the greater fraction of the total star formation occur- 
ring in non-starburst galaxies. This is supported by the 
results presented in Section ^31 which suggest that the 
decrease in psfr(z) with cosmic time is largely driven by 
the secular decay in SFR after earlier accretion events. 
That is, the current evolution in psfr{z) is largely driven 
by interior baryon physics rather than the merger of 
CDM dominated halos. This interpretation is consistent 
with the relatively small fraction of multiple ELG sys- 
tems and targets that look like recent mergers (15 out 
of 93 HiPASS pointi ngs included in SR I, see Paper I), 
as well as the work of lBell et al\ l|2005f) . In contrast, at 
z > 1 field galaxies tend to have a disturbed morphology 
suggesting a higher fraction of interacting and merging 
systems and very few regular disk galaxies can be dis- 
cerned. It is tempting to speculate that z ~ 1 represents 
the epoch where cosmic evolution transitions from being 
driven largely by CDM interactions to being driven by 
the internal self-regulation of star formation. 
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Fig. 1. — Fraction of total luminosity density per Hi mass decade. Dark symbols represent Ha, while lighter symbols denote R band 
luminosities. Circles represent the value for the individual SINGG galaxies. Diamonds and error bars are the average values and standard 
deviations of mean for each mass bin. All values are corrected for internal extinction. 
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Fig. 2. — Gas cycling timescale and histogram. Circles are values for individual galaxies. Diamonds and error bars are the average 
values and standard deviations of mean for each mass bin. The dashed line corresponds to the Hubble time (13.6 Gyr). Dark symbols are 
corrected for internal extinction, lighter symbols are uncorrected values for the same galaxies. 
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TABLE 1 
Local Ha cosmic star formation rate density measurements 



Survey 


Sources 


z 


log{PSFR(^)) 


log(pSFR(^)) 


l°g(PSFR(0)) 


log(psFR(0)) 


SINGG 


110 


0.01 


r, 1, +0.08 

-z.io _o.o9 


1 8n +0-13 

-l.»U _o og 


14 +008 


1 o-i +0.13 


UCM (1) 


264 


0.02 


-2.45 ± 0.20 


-1.92 ± 0.20 


-2.48 ± 0.20 


-1.95 ± 0.20 


UCM (2) 


79 


0.03 


-2.15 ± 0.11 


-1.60 ± 0.11 


-2.19 ± 0.11 


-1.64 ± 0.11 


SDSS (3) 


149660 


0.10 




-1.54 ± 0.07 




-1.66 ± 0.07 


FOCA (4) 


216 


0.15 


-2.21 ± 0.15 


-1.82 ± 0.06 


-2.39 ± 0.15 


-2.00 ± 0.06 


CFRS (5) 


110 


0.20 


-1.88 ± 0.06 


-1.49 ± 0.06 


-2.12 ± 0.06 


-1.73 ± 0.06 



Note. — Units for columns 4-7 are [Mq yr~^ Mpc~^]; Pspr(2) and psfb.{z) are 
the SFR density estimates without and with internal dust extinction corrections, 
respectively, evaluated at the mean redshift z of the survey. The corresponding 
z = rates are extrapo lated assuming psfr{z) = /6sfr(0) (1 + z)^ 



Gallego et aB Im 



Perez-Gonzalez et al. (20& 
Brinchmann et al. (2004,) 

Sullivan et al. ( 2000) 

Tresse fc Maddox ( 199^ 







TA 
Hi mass 


BLE 2 

FUNCTIONS 








HiMF 
(1) 




a 
(2) 


log(A4*) 
(3) 


e, X 10^ 

(4) 


log{psFR) 
(5) 


log(PHl) 

(6) 


This paper 


(2002) 


-1.41 ± 0.05 

-1.37 ± 0.03 

-1.30 ± 0.08 

-1.53 


9.92 ± 0.04 

9.86 ± 0.03 

9.85 ± 0.06 

9.94 


3.9 ± 0.7 

4.9 ± 0.7 

7.5 ± 1.7 

4.72 


1 sn +0.13 

-l.»U _0 Qg 

+0.15 
!■' ' -0.11 
1 62 +013 
-, CO +0.14 

-1.08 _p ;^(J 


7.71 ±0.03 

r, yo +0.06 


Zwaan et al. C2003) 
Rosenberp & Schneider 


7.88 ±0.02 
7.91 ±0.01 







Note. — Column descriptions [units]: (1) Source reference. (2) 
Schechter fit power- law constant. (3) Schechter fit characteristic HI mass 
[AI0]. (4) Schechter fit normalization [Mpc~^ dex~^]. (5) Star formation 
rate density [Mq yr"^ Mpc"^]. (6) Hi mass density [Mq Mpc"^]. 



TABLE 3 
Cosmic star formation as a function of mass 



log(A^Hi/X0) 


N 


tgas 


PSFr per log(XHi/A^©) 


In per log{Mm/MQ) 




(1) 


m 


(3) 


(4) 


(5) 


7.0- 


8.0 


4 


5.9 ± 0.5 


(1.68 ± 1.31)xl0-3 


(5.53 ± 3.06)xl036 


8.0- 


8.5 


13 


6.1 ± 0.4 


(2.88 ± 3.23) X 10-3 


(5.15 ± 2.69) X 1036 


8.5- 


9.0 


20 


21.4 ± 3.1 


(1.41 ± 0.61)xl0-3 


(5.66 ± 2.52) X 10315 


9.0- 


9.5 


18 


9.2 ± 2.3 


(5.96 ± 1.38)xl0-3 


(41.9 ± 33.1)xl036 


9.5- 


10.0 


21 


6.6 ± 1.1 


(10.2 ± 2.3)xl0-3 


(40.2 ± 8.2) X 1036 


10.0 


- 10.5 


14 


5.9 ± 1.0 


(7.28 ± 1.98)xl0-3 


(32.5 ± 5.2) X 1036 


10.5 


- 11.0 


2 


6.9 ± 0.7 


(0.71 ± 0.17)xl0-3 


(2.48 ± 0.13)xl036 



Note. — Column descriptions [units] : (1) Himass range. (2) Number of galaxies 
within Hi mass range. (3) Gas cycling timescale [Gyr]. (4) Star formation rate 
density contribution per decade of Hi mass [A^q yr"^ Mpc"^ dex~^]. (5) R band 
luminosity density contribution per decade of Hi mass [erg s~i A"'^ Mpc""^ dex~i]. 
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TABLE 4 
Error budget for luminosity densities 



Uncertainty 


Uncertainty in log 


(luminosity density) 




I'r 


'fl 


'h. 


'Ha 


Random Errors 










Hi mass function 


±0.029 


±0.027 


±0.034 


±0.029 


Sampling 


±0.062 


±0.073 


±0.071 


±0.061 


Sky subtraction 


±0.001 


±0.001 


±0.002 


±0.002 


Continuum subtraction 






±0.010 


±0.015 


Flux calibration 


±0.008 


±0.008 


±0.013 


±0.012 


Nil correction 






+0.004 
-0.006 


+0.005 
-0.008 


Internal dust extinction 




+0.053 
-0.013 




+0.107 
-0.016 


Total Random 


±0.069 


+0.094 
-0.080 


+0.078 
-0.083 


+0.127 
-0.075 


Systematic Errors 










Nil zeropoint 






±0.002 


±0.002 


Internal dust zeropoint 




±0.003 




±0.006 


Distance model 


±0.014 


±0.017 


±0.028 


±0.033 


Total Systematic 


±0.014 


±0.017 


±0.028 


±0.033 



TABLE 5 
Breakdown of lua and In by galaxy parameters 



Percentile 


10% 


25% 


50% 


75% 


90% 




(R/Ua) 


(R/Ra) 


(R/Ra) 


(R/Ua) 


(R/Ua) 


log{Mm [Mq]) 


8.19/7.96 


9.25/9.08 


9.52/9.58 


9.97/9.97 


10.13/10.12 


log(A4dyn [Mq]) 


9.55/9.27 


10.34/9.87 


11.07/10.67 


11.42/11.16 


11.45/11.42 


Mr [ABmag] 


-16.86/-16.63 


-19.64/-18.40 


-21.68/-20.77 


-22.49/-21.84 


-22.69/-22.42 


log(SFR [M@ yr^i]) 


-1.19/-1.05 


-0.27/-0.33 


0.19/0.17 


0.45/0.74 


0.89/1.04 


A (Ha) [ABmag] 


0.34/0.32 


0.66/0.51 


1.11/0.92 


1.38/1.18 


1.45/1.36 


\og{r,{R) [pc]) 


3.01/2.88 


3.47/3.37 


3.61/3.61 


3.73/3.73 


3.89/3.88 


log(re(Ha) [pc]) 


2.90/2.73 


3.42/3.17 


3.69/3.54 


3.84/3.83 


4.07/3.99 


log(re(Ha)/re(/?)) 


-0.19/-0.49 


-0.06/-0.11 


0.06/0.02 


0.14/0.09 


0.22/0.14 


r9o{R)/re(R) 


1.96/1.89 


2.06/1.99 


2.47/2.23 


3.26/2.91 


3.74/3.29 


log(EW5o(na) [A]) 


0.53/1.22 


1.05/1.43 


1.40/1.64 


1.64/1.87 


1.87/1.97 


l^eiR) [t] 


22.26/22.28 


21.42/21.44 


20.35/20.43 


19.35/19.97 


19.04/19.33 


log(Se(Ha) [t]) 


-2.57/-2.55 


-2.31/-2.28 


-2.19/-1.72 


-1.58/-1.28 


--1.02/-0.68 



Note. — 
t: units are [ABmag arcsec"^] 
*: units are [Mq yr^^ kpc^^] 



